Supermodes of photonic crystal CCWs and multimode bistable switchings with 

uniform thresholds 
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Photonic crystal (PC) coupled cavity waveguides (CCWs) with completely separated eigenfre- 
quencies (or supermodes) are investigated. Using a coupled mode theory, the properties of the 
supermodes, such as the electric field profiles, the eigenfrequencies (including the central frequencies 
and corresponding linewidths), and the quality factors (Qs) are expressed in very simple formulas, 
and these results agree well with those of "exact" numerical method. We also discuss the great 
differences between the supermodes and continuous modes, which are described by tight binding 
(TB) theory. Then, the properties of supermodes are used to investigate their potential applications 
in multichannel (covering the whole C-band) bistable switchings, and bistable loops are obtained 
numerically. We also predict and verify that all the thresholds are uniform inherently. 
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INTRODUCTION 

Photonic crystal (PC) 0,12 coupled cavity waveguides 
(CCWs) 13, Q, which are formed by series of coupled 
point defects in otherwise perfect PCs, have been inves- 
tigated intensively using tight binding (TB) theory in 
the last several years IM IE S IS 0' ^^'^ propaga- 
tion mechanism of hopping between neighboring cavities 
has been clearly understood 0. Generally speaking, a 
periodic CCW possesses a continuous transmission band 
of 17o(l - Aq:/2 -I- 1 - Aa/2 + \k\), where flo is the 
eigenfrequency of an individual cavity and n, Aa are 
two related overlap integrals |3j01- For convenience, the 
modes of the continuous band derived from TB theory 
are called TB-modes in this paper. MarpyDotential ap- 
plications, such as broadband dropping 11], broadband 
splittin g 12,, J_^] , broadband optical switching and limit- 
ing [l3.ll5|Tand ultrashort pulse transmission 16] are all 
based on the continuous transmission bands (TB-modes). 

In practice, however, the number of the cavity N is 
finite, and Born-von Karman periodic condition is used 
0, and N discrete modes are obtained, which satisfy 
the dispersion relation of 0, @] ■ 



w(fcm) = r^o [1 — Aa/2 + KCOs{kmR)] 
with fc„i ~ 2rmr/L, m — 0, ■ ■ ■ , N — 1 



(1) 



Where R is the distance between two neighboring cav- 
ities, and L is the total length of the system. When 
the coupling strength between the cavities are designed 
carefully, and some criteria are satisfied ^3] i ^ quasi-flat 
spectrum with N small dip s formed by the N discrete 
modes may be obtained [l7j]. Actually, the broadband 
operations mentioned above [ill . lla . Q, 0] are all 
realized using the quasi-flat transmission band. 

On the other hand, one can also design a A^-cavit 
CCW to make the N modes separated completely 
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FIG. 1: One dimensional photonic crystal CCW structure 
considered in this paper. The individual cavity of C is repre- 
sented by C = ababcbaba. Where a, b and c are three differ- 
ent kinds of medium layers with refractive indexes of Ua, Uf, 
and lie, and physical thicknesses of Ao/(4na), Ao/(4nb) and 
Ao/(2nc), respectively. The layer a is the same as layer a 
except that the thickness is a/2. Ao is a reference wavelength. 



(these separated modes are called supermodes in this pa- 
per and the reasons are given below), rather than to 
form a quasi-flat band. We find that the predictions 
of TB theory become inexact (as shown below) in this 
case. Although the supermodes of a A-cavity system 
have been noticed early ^3] ? the properties of them are 
not investigated intensively (such as the mode profiles 
and linewidths of each supermodes) , which prevents them 
from wide application. 

According to the criterion developed in Ref. 0], one 
dimensional (ID) CCW with discrete supermodes is de- 
signed, as shown in Fig. ^ Although some numeri- 
cal methods, such as the finite difference time domain 
method (FDTD) |l3], transfer matrix [UliOl can be used 
to extract the properties we expected, we prefer to un- 
derstand them from the viewpoint of physics and express 
them in analysis and simple formulas. 

In this paper, a general coupled mode theory is pre- 
sented to analysis the supermodes of A^-cavity CCW sys- 
tems, as shown in Fig. ^ The eigenfrequencies, includ- 
ing both the central frequencies and the corresponding 
linewidths (and also the quality factors), as well as the 
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mode profiles of the supermodes are obtained, and agree 
well with "exact" numerical results (obtained using stan- 
dard transfer matrix method). Subsequently, these prop- 
erties of supermodes are used directly to design a poten- 
tial application of multichannel bistable switchings in the 
whole C-band. Using the results of coupled mode theory, 
we prove that the thresholds of the multichannel switch- 
ings are low and uniform when Kerr media are carefully 
introduced into the cavities. 



COUPLED MODE THEORY ANALYSIS OF 
SUPERMODES 

In the coupled mode theory presented in this paper, the 
electric field (r) of the entire coupled cavity system is 
expressed as a linear superposition of the modes of the 
N cavities: 



N 



(2) 



where A„ {n = 1, • • • , N) are complex coefhcients that 
determine the relative phase and amplitude of the cav- 
ities. E^^; and Esi(r — nRez) are the fields of the entire 
coupled system and those of individual cavities centered 
at uRgz respectively, is the direction of the cavities 
being aligned, to and 17 are the allowed frequency of the 
coupled system (is unknown now) and the frequency of 
an individual cavity (is already known). We normalize 
Ef2(r) to be unity according to /e"(r)E(r) • E(r) = 1, 
with e'^(r) the dielectric function of an individual cavity. 
Eq. ((2Jl is very similar to the linear superposition field 
used in TB theory, which reads 



Eif(r) = ^B„En(r-ni?e^ 



(3) 



n=l 



where the superposition coefficients are i?„ = 
Eq eicp{—inK R) . However, one may find great differ- 
ences between them. In TB theory (Eq.(|31l), the superpo- 
sition coefficients are i?„ = Eoexp{—inKR), which have 
the same modulars of \Eo\ for all the cases of n from 
1 to N, and the relative phases between them are also 
determined. However, in Eq.Q, the coefficients of An 
are arbitrary complexes, and the amplitudes and phases 
of them may be greatly different for various of n. This 
is one of the most important differences between super- 
modes and TB-modes. 

Substituting Eq.Q into the simplified form of 
Maxwell's equations of: 



V X V X E = e(rV 



(4) 



Then, we operate both sides of the resulting equation 
from left using the operator of / ci'E(r — mRez)-, one 



can obtain a group of coupled equations: 

N 

Am ^ ^ ^ CninAn — 

where the coefficients Cmn are defined as 



(5) 



Crn.n. — 



_ (1 + Aa„)w2 
(}„in ~ j dre{r — nRez)E{r — mRe^) ■ E(r — nRe 



(6) 



(7) 



amn= J dre{r)'E{r - mRe^) ■'E{r - nRe^) (8) 

Aa,„ = y"drAeE(r - mfle^) • E(r - TOiJe^) (9) 

(with Aem(r) — e(r) — e(r — mRr)) 

Eq. (jSJ is very similar to the governing equations of 
coupled waveguide arrays [2ll |. or phase- locked injection 
laser arrays [23 . E3 | , where, the solutions are called 
supermodes of the waveguide (or laser) arrays. Similarly, 
we name the solutions of Eq.lSJ as the supermodes of 
the CCW. Generally, Eq. (O can not be solved in closed 
form, however, it can be simplified and solved for some 
very special cases. For example, when only the nearest 
neighboring coupling are considered and the cavities are 
uniformly spaced. Then the coefficients are simplified to 



Cn 



1^2 - (1 + Aa)w2 



C, n — m ± 1; 
0, others. 



(10) 



Here, we have used the relations of Pm,m±i — 

am,m±i = a and Aa„ = Aa. 

Then, using the same method as used in Ref.|2ll l22l 
E3. 01 and proper boundary conditions, one obtain the 
solutions of the Eq. 0. Generally, there're N solutions 
(supermodes) for a A^-cavity CCW system. For the Lth 
supermode, the linear superposition coefficients A^ and 
C^, respectively, are: 



A^ sm{ne^) 
1 

~2cos(6l^) 
Lit 



1, 



N + 1' 



1, 



(11) 
(12) 

(13) 



Where A^ is a constant and is determined from the nor- 
malization condition of 



dre{r) 



N 



.ra=l 



JV 



^ A„Ef^(r- ne^-) 



.n=l 



(14) 

After a simple algebra process and using the normaliza- 
tion condition of the individual cavity modes, one can 
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FIG. 2: Comparison of the eigenfrequencies derived from the 
coupled mode theory and from numerical method for the 
structure of Fig. Q The refractive indexes are Ua ~ 2, 
rib = 1 and ric — 2 respectively. Number of cavities is = 16. 
(a) Transmission spectrum calculated using one dimensional 
transfer matrix method, (b) Comparison of the central fre- 
quencies of numerical results (Fig. EI a)) with those derived 
from Eq. lO with parameters of a = 0.00499, /? = -0.01497 
and Aa = (See text for details), (c) Comparison the 
linewidths and quality factors of the supermodes derived from 
transfer matrix method (Fig. |Ha)) and coupled mode theory 
(Eq. linj) with the parameter of Qo = 2.6347 x 10*. 



obtain the constant of is = + 1), which is 

independent of the mode number of L. 

According to Eq. JBJ and Eq. H13() , one can obtain the 
iV eigenfrequencies of the supermodes: 



CO 



C^Ac 



L = 1,2, 



,N (15) 



Eq. lO and Eq. (|15|l are the main results of the 
coupled mode theory discussed above. We want to 
point out that in Eq. (|12|l . tends to infinite when 
0^ — mn + 7r/2 (m is an integer). For example, when 
N = 5 and L = 3. However, this does not mean that the 
coupled mode theory is invalid in this case. Because the 
parameter of is a function of frequency of the super- 
mode, a and /3 (see Eq. (jSJ), but not an physical quality. 
The values of physical qualities of A^ (Ea. l(TT|l ) and uj^ 
(Eq. H15|l ) are both finite and correct when compared 
with numerical results. 



Actually, the method shown above is not new, and it 
has been considered in Ref. for the special cases of 
N — 2 and 3. When iV = 2, one can derive from Eq. l(T^ 
that C^'^ = =f1, and then from Eq. pSI) . one can obtain 
the frequencies of the supermodes and the superposition 
coefficients: 



l±f3 



1 + Aa±a 



1 

72' 



1 



V2 
(16) 

Clearly, the superposition coefficients and the corre- 
sponding frequencies in Eq. Hl()|l are the same as the 
results of Eq. (2) in Ref. (J). Similarly, one also can 
calculate the frequencies and fields for the case of = 3, 
which are also the same as the results in Ref. 0] . (In Eq. 
(3) of Ref. 0], the Aa is regarded as negligible small). 

In order to show the power of the coupled mode the- 
ory discussed above, we consider a one dimensional CCW 
with a large number of cavities (i.e., N = 16), as shown 
in Fig. |T] Using the standard transfer matrix method 
(TMM) [l9|, one can easily obtain the transmission spec- 
trum, as shown in Fig. HJa). Using Eq. H15() . H12() and 
one also can calculate the 16 eigenfrequencies, and 
the results are shown in Fig. [2Jb). Clearly, the theoret- 
ical results agree well with the numerical results. Here, 
the parameters are set to be a = 0.00499, /3 = -0.01497 
and Aa = 0, which are obtained according to the follow- 
ing two steps: Firstly, we numerically obtain the two 
eigenfrequencies of (w^'V^^o) = 0.990017, 1.009983 for 
the case of A = 2. Secondly using Eq. (|16|l . we calcu- 
late the above parameters of a, f3. Here we assumed that 
Aa is negligible compared to a and /3. On the contrary, 
when the values of a, /3, Aa and km are substituted into 
Eq. much larger errors of uj (km) from the numerical 
results are found. Therefore, the coupled mode theory 
gives more exact results than the TB theory in the case 
of discrete supermodes. 

The quality factors (Q's) of the supermodes can also 
be derived from the results of coupled mode analysis. 
The Q's of the Lth supermode is defined as — 
Lu^W^/AW^, where = ^ Jy drE • D is the total en- 
ergy stored in the coupled cavity system of the Lth super- 
mode, and AW^ = -dW^/dt is the energy dissipated. 
According to the normalization condition of Eq. H14|l , one 
can easily find that the total energy of all the A^ super- 
modes are the same, i.e, = Wq. For the CCWs com- 
posed of lossless media, the dissipation of energy of the A^ 
supermodes are all due to the coupling out of the system 
through the 1st and A^th cavities. Therefore, the energy 
dissipation of AW^ is proportional to the energy stored 
in the first and last cavities, i.e., AW^^ = AW^o sin^(6'^). 
Here, we have used the results of Eq. Hll() . Then, one 
can find the quality factor of the Lth supermode: 



AW^ 



UJ 



Wo 



(17) 
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Here, Qo = ^oWq/AW, and we have used the the fact 
of u!^ w r^Q. The results of Eq. (|17|l have been shown in 
Fig. ^c) with the parameter Qq = 2.6347 x 10''. 

On the other hand, we can also derive the by find- 
ing the values of = uj^ / Aujj^f^^ from the numerical 
results (Fig. Ill^a)), where Aw^^^,^ is the full width at 
half maximum of the Lth supermodes. The results are 
also shown in Fig. IHc), and they agree well with the 
results of Eq. lfT7| . 

Fig. 121 shows the field profiles of the first half super- 
modes of u!^''"'^. For the case of clarity, the modes of 
^9,- -,i6^ which are very similar to the modes of w^'^' "'^ 
respectively, are not shown. Clearly, the results of cou- 
pled mode theory of Eq. are agree well with those 
of the transfer matrix method. The amplitudes of indi- 
vidual cavities of the Lth supermodes lie on an envelope 
function of sin(-^q^). Therefore, the amplitude of some 
cavity modes may be tend to zero, while others reach to 
maxima. In the TB model, however, the amplitude of all 
the cavity modes are the same. 

From the discussion above, one can find that the cou- 
pled mode theory presented above describes the main 
characteristics of the supermodes very well, provided that 
the parameters of a, (3 and Qq are given. The super- 
mode states of a CCW are very different with the states 
of TB-modes, which are described by the tight binding 
theory 3, 4]. For the TB-modes, the transmission band 
is determined by Eq. (Q, while in the supermodes, the 
central frequencies are determined by Eq. H15() . In the 
TB modes, the localizations in each cavities are the same, 
while in the supermodes, the localizations change greatly 
according to a simple sine function. In the next section, 
one of the possible applications of supermodes of multi- 
channel bistable switchings is proposed and analyzed us- 
ing the results of coupled mode theory given above. 
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FIG. 3: (Color online) The electric field profiles of supermodes 
calculated using transfer matrix method (blue lines) and the 
envelope of superposition coeflicients (red lines with circles) 
shown in Eq. IIH . The distribution of the refractive index is 
also shown for clarity. 
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FIG. 4: (Color online) (a) Bistable switching loops and the 
thresholds of the 16-channel switchings, (a) Switching loops 
of i^l^ and those of i^f^ '^^ are not shown for clarity, 
which are similar to Normalized unit of light inten- 

sity of 712/ is used with n2 the Kerr nonlinear coefficient and 
/ the light intensity, (b) the thresholds of the 16-channel 
switchings. 



MULTI-CHANNEL BISTABLE SWITCHINGS 
WITH UNIFORM THRESHOLDS 

One of the potential applications of the supermodes is 
multichannel optical bistable (OB) switching, which is a 
key component in all optical information systems. The 
OB switching have been widely investi gat ed in photonic 
crystals with Kerr defects El IM El IH , of which 
the refractive index changes with local light intensity, 
i.e., n = uq + 77.2/ with uq the linear refractive, n2 the 
nonlinear Kerr coefficient and / the local light intensity. 
However, most of the researches are focused on a single 
frequency operation . In a wavelength divi- 
sion multiplexer (WDM) system, there're always need a 
multi-channel switching for all the working channels. We 
find that the supermodes of CCWs are very suitable for 
this function [28j. As an example, we investigate the 
systems of 16-cavity CCW structures shown in Fig. 

According to the OB switching theory, the shift of the 
eigenfrequencies with the changes of the dielectric con- 
stants of the cavities is a key factor for OB operation. 
Using the perturbation theory , one can find the shifts 
of eigenfrequencies of supermodes with the change of di- 
electric constant Ae(r): 

'.■^/A6(r)|E^.(r)p 



At^^ = — 



where E = ^«E(r 
permodes of the unperturbed CCW 



2 /e(r)|E^.(r)P ^^^^ 
nez) is the electric fields of su- 
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When only the 1st and A''th cavities are perturbed by 
the same amount of Ae, then using the relation of Eq. 
((TT|l . ((TH) and one can derive that 

Ac^^ = -^ / A6|E„.(r)|2dr = Ac^o|sin(0^)|2 

^ JCi,Cn 

(19) 

Where, Acjo = -[cj^Ae/(iV + 1)] d-jE^^ (r) + E^jl (r - 
iVi?ez)p. On the other hand, according to Eq. (|17|) . we 
find that the FWHMs of the supermodes are w^^^,^ = 
(w^/Qo) sin^(0^). Therefore, the shifts of the frequen- 
cies are proportional to their FWHMs respectively, i.e., 
Ao;-^ — {Au)oQo/^o)i^fwhm- When the incident frequen- 
cies Wj^j of the bistable switchings are tuned from the 
central frequencies according to w^^^ = uj^ — yujj^f^^, the 
thresholds of the multi-channel bistable switching are ex- 
pected to be uniform. 

Using the nonlinear transfer matrix method |30l |. and 
setting y = 2, we obtain the bistable switching loops, and 
the results are shown in Fig. 0] In Fig. 01 we have used 
the normalized intensity of J^]. We can see that the 
thresholds of the switchings are almost the same. Except 
for the 1st and 16th channel (In practice, the 1st and 
TVth channels may be tuned slightly in order to switch 
them with approximately the same thresholds as other 
channels), the maximum of the thresholds is 4.8 x 10~^ 
(the channel of uij^^) and the minimum is 4.5 x 10~^ (the 
channel of ajf„), and the relative difference is about 6%. 

When the thicknesses of a, b layers are 193. 8nm and 
387. 5nm respectively, the central frequency of Aq is at 
about 1.55/Ltm, and the multi-channel frequencies cover 
the C-band of optical fiber communication entirely. For 
a Kerr nonlinear coefficient of n2 = 2 x 10~^^cm^ /W (a 
value achievable in many nearly instantaneous nonlinear 
materials), the thresholds are about 4.6 x \Q~^W/cm?, 
which is much smaller than those of the switchings stud- 
ied before 

CONCLUSION 

In summary, we have analyzed the supermodes of 
CCW systems, which are different from the TB-modes 
and are also important operation states of CCWs. Using 
the coupled mode theory, the eigenfrequencies, including 
the centers and the FWHMs, quality factors, and mode 
profiles of the supermodes are formulated in very simple 
forms. And they agree well with exact numerical results. 
We also discussed the great difference of the supermodes 
with the quasiflat TB-modes. We investigated one of 
the potential applications of the supermodes, which is 
a 16-channel bistable switching with uniform thresholds 
covering the C-band of optical fiber communication. The 
results show that the thresholds are low and uniform. 
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